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We show that the vast majority of all pure states featuring a common expectation value of some 
generic observable at a given time will yield very similar expectation values of the same observ- 
able at any later time. This is meant to apply to Schrodinger type dynamics in high dimensional 
Hilbert spaces. As a consequence individual dynamics of expectation values are then typically well 
described by the ensemble average. Our approach is based on the Hilbert space average method. 
We support the analytical investigations with numerics obtained by exact diagonalization of the 
full time-dependent Schrodinger equation for some pertinent, abstract Hamiltonian model. Further- 
more, we discuss the implications on the applicability of projection operator methods with respect 
to initial states, as well as on irreversibility in general. 
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In its broadest sense the term typicality may be ex- 
plained as follows [J, If a set of states specified by 
some common feature (e.g., drawn according to the same 
distribution, sharing the same energy, etc.) yields a very 
narrow distribution of some other feature (e.g. some ob- 
servable, etc.), then there is typicality. The concept of 
typicality as key to the occurrence of standard statistical 
equilibrium behavior (as opposed to ergodicity, mixing, 
etc.) especially in quantum mechanics has recently been 
established in various works [1, |3j S 0] ■ 

One important implementation has been presented in 
There it is shown that a majority of pure states of 
a compound system from some narrow energy shell yield 
almost the same reduced density matrix for a small sub- 
system. This happens to be the canonical equilibrium 
state in case of weak interactions and generic environ- 
ment spectra [3]. 

In another typicality-based investigation it is 
demonstrated that states drawn according to a certain 
type of probability distribution in Hilbert space fea- 
ture very similar quantum expectation values (QEV's) 
of generic observables. Even more detailed results exist 
for the special case of a uniform distribution of normal- 
ized, pure states, i.e., a distribution which is invariarit 
under all unitary transformations in Hilbert space 
Using the Hilbert space average method (HAM) 
one finds that the "Hilbert space average" (HA), i.e., the 
average of the QEV's of an observable D w.r.t. the above 
distribution, is given by 



HAK^Ii^lV-)] = = ci 



(1) 



and the corresponding "Hilbert space variance" (HV) by 
YWmom := HA[((^|i?|V') -HAKV-li^l^)])'] 

n + 1 \ n \ J I n+1 

cf . . Here n denotes the dimension of the correspond- 
ing Hilbert space and ci := TrjD'l/n describes the i-th 



moment of the spectrum of D. Thus C2 — is the spectral 
variance of D. Throughout this paper we focus on observ- 
ables the low spectral moments of which do not change 
(significantly) under physically reasonable "upscaling" . 
Pertinent examples are, e.g., a component of a specific 
spin in a system which is upscaled by adding more and 
more (interacting) spins, the occupation number of some 
momentum mode in an interacting many-particle system 
which is upscaled to comprise more and more momen- 
tum modes, or in general any local variable embedded in 
a growing system. For any such observable one may con- 
clude from (21) that the Hilbert space variance, i.e., the 
width of the above distribution of QEV's vanishes with 
growing dimension n. In this sense bound observables in 
large systems yield typical QEV's. 

In this paper we turn towards the typicality of dynam- 
ics of QEV's. In short, we demonstrate in the paper 
at hand that pure states from a set {!(/>)} featuring a 
common QEV of some observable A at some time t, i.e. 
{(t)\A{t)\4)) = a, most likely yield very similar QEV's at 
any later time, i.e. {(l)\A{t + t)|(^) « {4i'\A{t + t)|(/)') 
(with 10), |(/))' both being states from the above set). We 
present some analytical derivations based on the HAM, 
in particular on Eqs. ([T]and[2]) and we additionally sup- 
port the results with numerical calculations. Finally, we 
discuss what consequences arise for the validity of pro- 
jection operator methods (Nakajima-Zwanzig (NZ), etc. 



(lOl . Illl |12| |) w.r.t. initial states and the corresponding 
inhomogeneities. Furthermore, we comment on the ir- 
reversibility of QEV's corresponding to individual pure 
states. 

We specify our considered observable A only by the 
above mentioned moments, c^, and specialize without 
substantial loss of generality to observables which are 
trace- free, ci — 0, and normalized to C2 = 1. Further- 
more we require the ci with z = 2, 8 to be of the order 
1. Next, we introduce an ensemble of pure states 
which is characterized as follows: All its states must fea- 
ture the same QEV of the observable A, {(l)\A\(j)) = a, 
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must be normalized {{4>\4>) = l)j f^^id uniformly dis- 
tributed otherwise. That means the ensemble has to stay 
invariant under all unitary transformations in Hilbert 
space that leave the expectation value of A unchanged, 
i.e. those transformations that commute with A, or, con- 
cretely, transformations of the form e'^, with [-8, A] = 0. 
This specifies the most general ensemble consistent with 
the restriction that all its states should yield a given a. 

For the following calculations we further introduce 
some kind of "substitute" ensemble which is much 

easier to handle. As will be shown below, this ensemble 
approximates the exact ensemble {\4>)} described above 
very well for large Hilbert spaces. 

The ensemble {|w)} is generated by 



(3) 



where are pure states drawn from a uniform distribu- 
tion of normalized states without further restriction as 
described above ([!]). d is some small parameter which 
describes the deviation from the "equilibrium" ensemble 
{IV')}- Since it is essentially the operator A itself that 
generates from the entirely uniform distribution, 

{|w)} is invariant under the above uniform transforma- 
tions that leave a invariant. 

The construction ([3|) allows for an evaluation of mo- 
ments of the distribution of {ijj\C\bj) based on results on 
moments of the distribution of {ip\D\^)^ or concretely 

nk[{uj\C\ujY] =Y^k[{'<p\D\ijY] 

with D = —^{l + dA)C{\ + dA). (4) 
1 + (1^ 

(Of course the average on the l.h.s. corresponds to the 
substitute ensemble {|cj)} while the average on the r.h.s 
is based on the completely uniform ensemble {IV')})- Ex- 
ploiting this, average and variance of (aj|C|w) may be 
evaluated with the help of (|ll2p . 

To assure that the ensemble {l^i-")} indeed approximates 
the ensemble {If/))}, in the limit of large n, we evaluate 
the following four quantities 



HA[(a;|a;)], 
HV[(a;|cc.)], 



Y{k[{uj\A\u)] , 
YW[{Lj\A{t)\uj) 



(5) 



where A{t) denotes the time dependence according to the 
Heisenberg picture. (For clarity: the results are given in 
Eqs. ®, 0, ® and Cni).) 

The states are not exactly normalized which would 
render them unphysical, of course. However, one finds 
from ^ and ^ (by implementing C ~ \) that 

HA[(w|w)] = 1 . (6) 

By exploiting ([2]) and ([4]) one finds analogously for the 
variance 



HV[(w|w)] 



1 4rf2+4d3c3 + rf''(c4-l) 



n+1 



(7) 



As defined above, the q are of the order 1, i.e. the HV 
of the norms scales with 1 /n and becomes small for large 
Hilbert spaces. Therefore, the vast majority of the states 
\uj) are approximately normalized for large n. 

The average of the QEV's of A w.r.t. the ensemble 
{|a;)} (which is meant to correspond to the above a) is 
calculated by exploiting ([T]) and (|4|) (by implementing 
C^A) 



Yik[{u\A\uo) 



2d + d^cg 



(8) 



That is, the mean QEV can be adjusted through the 
choice of the parameter d. However, the replacement 
ensemble is restricted on expectation values not too far 
away from zero (i.e. the average expectation value of 
the "equilibrium" ensemble {IV')}) because by sweeping 
through all possible d not all possible expectation values 
up to the maximum eigenvalue of A are reachable. 

The evaluation of HV[(Li;|74(t)|a;)] turns out to be some- 
what more complicated, since we, in general, cannot fully 
diagonalize the Hamiltonian and thus do not know A{t) 
in detail. However, we are able to perform an estimation 
for an upper bound. For this purpose we make use of 
the Hilbert Schmidt scalar product for complex matrices 
defined as [X^Y] := Ti{X''Y}. Thus, one can formulate 
a Cauchy-Schwarz inequality of the form 



Tt{X'<Y} < A/Tr{XtX}Tr{rty} 



(9) 



Particularly, one obtains Tr{A{t)A} < Tt{A'^}. Evaluat- 
ing HV[(w|A(t)|w)] based on ^ and g]) (by implement- 
ing C = A{t)), realizing that Tr{_D}-^ is always positive 
and repeatedly applying ([9]) yields the inequality 



RV[{Lo\A{t)\oj)] < ■ 
1 -I- + Gd'^Ci + M^y/c^{/c^ + d^yjc^ 

(1+^2)2 



(10) 



Again, since the Ci are of the order 1, the upper bound 
deceases with 1/n. Thus, the variance (fTO)) becomes 
small for large Hilbert spaces, just like the variance of 
the norms ([7]). This result yields two major direct impli- 
cations. 

First, if one evaluates (|10p at f = 0, one finds that the 
majority of the states \lo) feature approximately the same 
QEV of the observable A for large n. From this property 
together with the result that the states are nearly 
normalized one concludes that the replacement ensemble 
{jw)} indeed approximates the exact ensemble {|(/))} very 
well for large Hilbert spaces (with a ~ HA[(ix'|A|w)] as 
given in (O). 

Second, the upper bound from ((T0|) is valid for any 
time t. Thus, for large enough systems, the dynami- 
cal curves for a^[t) := (ci;|A(f)|w) of the vast majority 
of pure states from the initial ensemble {|w)} are very 
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close to each other and thus to the evolving ensemble 
average at any time t. Due to the similarity of {\uj)} 
and {|0)} this should also hold true for the "exact" en- 
semble {10)}- Thus, there is a typical evolution for the 
expectation values {<j)\A{t)\(j)) or, to rephrase, there is 
"dynamical typicality". This statement represents the 
main result of this paper. Particularly, this typicality is 
independent of the concrete form of the dynamics, which 
may be a standard exponential decay into equilibrium or 
something completely different. 

In the following we visualize these predictions for a 
model quantum system described by a Hamiltonian of 
the form H = Hq + V, where Hq is some unperturbed 
Hamiltonian with equidistant energy eigenvalues (num- 
ber of states: n = 6000, level spacing: AE = 8.33•10-^ h 
set to 1), and V some possibly but not necessarily small 
interaction. A is chosen as diagonal in the eigenbasis 
of Hq with equally many, randomly placed elements 1 
and —1. This is in accordance with the already men- 
tioned conditions on the moments Ci (ci = , C2 = 1). 
However, this specific form of the Hamiltonian and the 
observable A is not crucial for the main results of this 
paper concerning the typicality of expectation values. It 
is just an example for the numerical illustrations. We 
further calculate the dynamical curves of a^^t) as gen- 
erated by three concrete interactions V that are chosen 
to represent different generic types of dynamical behav- 
ior. Thus, the matrix elements of V in the eigenbasis 
of Hq (and A) are taken as (i.) random complex Gaus- 
sian numbers with Vij — 0, |Vijp = 2.25 • 10~* (small 
perturbation), (ii.) random complex Gaussian numbers 
with Vij = 0, — 6.25 • 10"^ (strong perturbation). 



(iii.) all identical constants, i.e., V^, = 2.25 • 10 



We 



obtain the dynamics of the ai^{t) by numerically solving 
the full time-dependent Schrodinger equation). These 
interactions give rise to three archetypical evolutions (i.) 
an exponential decay into equilibrium (Fig. 1(a)), (ii.) 



non-exponential decay into equilibrium (Fig. 1(b)), (iii.) 
no decay to equilibrium at all, even in the limit of many 
states and weak perturbations (Fig. 1(c)). The precise 
reasons for the emergence of these dynamics are beyond 
the scope of this text (for more details see [13|). For a 
clarification of the term " relaxation" in this context refer 
to 



1^. 



In each of the figures the dynamics of Qui (t) for three 
different pure states of the ensemble {|a;)} as initial states 
are displayed. Furthermore, the evolution of the ensem- 
ble average is shown (averaged over 100 states). To nu- 
merically generate {Iij-")}, the real and imaginary parts of 
the amplitudes of the states \tp) (cf. ^) w.r.t. the eigen- 
basis of Hq are drawn as independent Gaussian num- 
bers, similar to the Gaussian Adjusted Projected ensem- 
ble (GAP) (see [H, [H). d is chosen as 0.1. Essen- 
tially, one finds that, for all three interactions, the dy- 
namical curves of the pure states are close to the average 
curve for all displayed times t. Furthermore, the insets 



in Figs. l(a)|l(b)|l(c) show the numerical variance cr^(t) 
of the a^{t) (calculated for 100 states), which turns out 
to be always smaller than the upper bound given in (|10p . 
Thus, all these numerics illustrate and back up our pre- 
viously derived analytical results. 




FIG. 1: Dynamics of some expectation values correspond- 
ing to a set of initial states {|<a')}, i.e., {uj\A{t)\uj) = auj{t). 
The set is characterized by an "common" initial expectation 
value, i.e., all yield {u>\A{0)\lj} ~ 0.2. The figures illus- 
trate that the individual evolutions typically stay close to the 
average over the set (solid line). The insets show evolutions 
of the variances a^{t) which stay accordingly small and re- 
main below some analytically predicted upper bound (solid 
line), cf. (llOf) . The subfigures correspond to different Hamil- 
tonians, generating diflterent archetypical types of dynamics: 
a.) exponential relaxation, b.) non-exponential relaxation, 
c.) non-relaxing. 



We now address further implications. The mean QEV, 
i.e., essentially a, can alternatively be reformulated using 
the notion of a density matrix as usually done in the 
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framework of projection operator formalisms 

a = llA[{uj\A\oj)] = RA[Tr{A\uj){Lu\}] ^ Ti{ARA[\lu) {uj\]} 

(11) 

The HA[|Li;) {uj\] takes the role of the density matrix. Fur- 
ther evaluation gives (using the "substitute" ensemble 
{\co)}) (seeii) 



2 a2 



l + 2dA + d^A 
n(l + 



(12) 



For ensembles close to equilibrium, i.e., small d, which is 
fulfilled in the examples presented here, one can neglect 
the terms which grow quadratically in d. In this case, the 
density matrix takes approximately the same form as the 
initial state which is often used in projection operator 
calculations which aim at determining the dynamics of 
expectation values like a(t) ([T§]). There, for reasons 
given below, the (mixed) initial state is simply taken to 
be p(0) — l/n + cA such that c = a(0). That means, 
correct dynamical results from the projection operator 
methods based on the above initial state describe the 
dynamics of the ensemble average of {Iijj)}- 

From this point of view some consequences on the 
applicability of projection operator theories (NZ, time- 
convolutionless, Mori formalism etc.), which are stan- 
dard tools for the description of reduced dynamics, arise. 
These methods have in common the occurrence of an in- 
homogeneity in the central equations of motion that typ- 
ically has to be neglected in order to solve them. Gener- 
ally, the inhomogeneity depends on the true initial state, 
it, however, vanishes if the true initial state indeed is of 
some specific form determined by the pertinent projec- 
17|. For the above mentioned case the above 
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p(0) is exactly of that form, which means the dynamics 
of the ensemble are equal to the dynamics generated by 
the pertinent projected equation of motion without the 
inhomogeneity. However, the evolution of the ensemble is 
typical, this implies that the inhomogeneity, as generated 
by most of the true initial states, should be negligible. 

On the other hand, there are investigations in the field 
of open quantum systems, e.g., 18 1 and [l^, suggesting 



that the true initial states may have an utterly crucial in- 
fluence on the dynamics, such that, e.g., some correlated 
initial states may yield projected dynamics which are en- 
tirely different from the ones obtained by corresponding 
product states. 

Nevertheless, to rephrase, the results of this paper indi- 
cate that in the limit of large (high dimensional) systems 
the inhomogeneity should become more and more irrel- 
evant in the sense that the statistical weight of initial 
states, which yield an inhomogeneity that substantially 
changes the solution of the projected equation of motion, 
should decrease to zero. Note that this does not contra- 
dict the concrete results of [l3| and [3]. 

The above results also shed some light on the rela- 
tion of the apparently irreversible dynamics of QEV's to 



the, in some sense, reversible dynamics of the underlying 
Schrodinger equation. If a mean QEV as generated by 
some initial non-equilibrium ensemble (pertinent density 
matrix) relaxes to equilibrium [3] (which can often be 
reliably shown [llj) , then for the majority of the indi- 
vidual states that form the ensemble, the corresponding 
individual QEV's will relax to equilibrium in the same 
way. Thus, for the relaxation of the QEV's, the ques- 
tion whether or not the initial ensemble truly exists is 
largely irrelevant. Of course, there may be individual 
initial states giving rise to QEV evolutions that do not 
(directly) relax to equilibrium, but, to repeat, for high 
dimensional systems, their statistical weight is low. 
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